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BAROCLINIC INSTABILITYWITH VARIABLE GRAVITY 
- A PERTURBATION ANALYSI S 

I. INTRODUCTION 


When the Space Shuttle becomes operational in the 1980's, some 
of the flights will take into orbit a laboratory known as Spacelab. Space- 
lab will have the environment of a laboratory on the Earth’s surface, 
except that when in orbit the value of effective gravity will be very 
small. A principal use of Spacelab will be to exploit this low gravity 
environment for science and technology. In particular, it will be 
possible to perform geophysical fluid flow model experiments in true 
spherical geometry. A radial dielectric body force, which is analogous 
to gravity, can be achieved over a volume of liquid held between two 
spheres. Such experiments cannot be performed in a laboratory on the 
Earth's surface because the dielectric body force cannot be made large 
enough to dominate the effect of terrestri^ gravity. A spherical model 
of the barocUnic, synoptic-scale processes of the Earth's atmospheric 
circulation is being designed for Spacelab. 

Since the last century, laboratory investigations seeking insight 
into large-scale circulations of the Earth's atmosphere have been per- 
formed. In these experiments a fluid was subjected to a horizontal 
temperature gradient and rotation; the temperature gradient modeled the 
polar-equatorial gradient and the rotation of the fluid container modeled 
the rotation of the Earth. The early experiments were qualitative in 
nature; systematic and quantitative developments date from the 1950's. 
Experiments in rotating cylindrical containers with heating at the rim and 
cooling at the base were carried out by Fultz et al. (1959). These 
experiments established the existence of two flow regimes, an axisymmetric 
flow for high heating and low rotation and an unsteady, non- axisymmetric 
flow for low heating and high rotation. Experiments in a rotating 
cylindrical annulus with differential radial heating were performed by 
Hide (1958), who discovered a steady, non- axisymmetric flow regime in 
which a jet flow meandered in a regular manner between the boundary 
cylinders. These annulus wave flows proved more amenable to systematic 
study than the cylinder flows (Fultz et al. , 1959; FowHs and Hide, 1965; 
Kaiser, 1970). 

Scaling of the governing equations for the large-scale features 
for both the atmosphere and the laboratory models leads to a set of 
approximate equations for both systems (Holton, 1972). These equations 
are known as the quasi- geostroptiic equations. Thus, the experiments 
are models of the synoptic-scale atmospheric flow. It is generally agreed 
that baroclinic instability is the process responsible for the large-scale, 
nonaxisymmetrlc flow in both the atmosphere and the models (Lorenz , 
1962). 



Although much has been learned from the cylindrical annulus 
experiments, further progress is inhibited by the geometrical difference 
between the cylindrical annulus and the Earth's thin, spherical- shell 
atmosphere. The obstacle to the extension of the cylindrical experiments 
to spherical geometry has been the difficulty of simulating a dominant 
radial body force. For a dielectric liquid contained between two con- 
centric spheres and subject to a voltage difference across the spheres 
in the manner of a spherical capacitor, it has been shown by Hart (1976) 
that a radial dielectric body force of the following form exists: 
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=oY (Ro'V 


( 1 ) 


where R. and are the radii of the inner and outer spheres , respec- 
tively; and are the dielectric constant and the density, respec- 
tively, at some reference temperature; n and y the temperature 
coefficients of dielectric constant and density, respectively; V is the 
voltage difference; and r is any radius between R. and R^. Values for 

the dimensions of the proposed Spacelab experiment and for the proper- 

2 

ties of suitable dielectric oils yield a typical value of g„ = 200 cm /sec 

£!j 

(Hart, 1976; FowHs and Fichtl, 1977). Thus, the need to perform these 
spherical model experiments onboard Spacelab becomes clear. Only in a 
low gravity environment will the dielectric body force dominate. 

The fact that effective Spacelab gravity will not be exactly zero 
is due to orientation rocket firings, crew movements, the vibration of 
Shuttle and Spacelab equipment , gravity gradient oscillations of the 
Shuttle, the deceleration in orbit due to atmospheric drag, etc. The 
combined effect of these sources will depend on many variables, but a 

2 -4 

typical value is estimated at 0.1 cm /sec , which gives a ratio of 5 x 10 
when compared to the value of the dielectric body force previously given. 
For special "quiescent" orbital flight periods, this ratio could be several 
orders of magnitude smaller. 

Design calculations are in progress for a spherical, baroclinically 
unstable model for Spacelab. Ultimately, a numerical, nonlinear, spherical 
model will be constructed, but to achieve understanding many theoretical 
calculations of increasing complexity have been undertaken. Since some 
of these calculations are also relevant to a wide range of research 
interests in theoretical dynamical meteorology and oceanography , they have 
been presented in publications (Fowlis and Arias, 1978; Geisler and 
Fowlis, 1979). 
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Equation (1) shows that the dielectric body force is a strong 
function of the radius. It was decided to calculate the effect of this 
variation on baroclinic instability. In Section II a quasi- geostrop hie, 
baroclinic instability problem in which gravity is a function of height is 
formulated. In Section III a solution is obtained using a perturbation 
procedure. In Section IV the results are discussed and the relevant 
information for the design of the Spacelab experiment is presented. The 
solutions are also relevant to other geophysical fluid flows in which 
gravity is constant but the static stability is a function of height. 


1 1. FORMULATION OF THE PROBLEM 


The quasi- geostrophic equations used to formulate the stability 
problem are well known and will not be re-derived here. Systematic 
derivations and discussions have been given by Pedlosky (1964) and 
Holton (1972). We consider the fluid to be contained in a rotating 
rectilinear channel with rigid horizontal and vertical boundaries. The 
horizontal radial and zonal directions are denoted by x and y, respec- 
tively, and z denotes the vertical direction. The nondimensional, 
linearized, incompressible, quasi- geostrophic equation for the perturba- 
tion pressure, p(x,y,z,t) is 



+ u 




9 q 9p _ 
9y 9 


( 2 ) 


where u(y,z) is the basic state zonal flow velocity in the x-direction and 
9q/9y is the meridional gradient of the potential vorticity of the basic 
state defined by 


ia = 

9y 



(3) 


F(z) is the rotational Froude number (inverse Burger number) defined by 
2 2 2 2 

F = f L /N H , where f is twice the rotation rate, and L and H are the 
respective distances between the boundaries in the y and z directions. 

1/2 

N(z) is the Brunt-Vaisala frequency defined by N = (yg 9T/9z) , 

where g is the acceleration of gravity and 9T(y ,z) /9z is the basic state 
vertical temperature gradient. The quantity 6 represents the scaled 
latitudinal derivative of f. 
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The nondimensional boundary conditions, which state that the 
normal component of flow at the boundaries is zero, are 

/ 9 9 N 9p 9u 9p _ n _ n 1 /•/i\ 

(■^ ^ Jx ) dz dTJT - > on z - 0, 1, (4) 

p = 0,ony = 0,l. (5) 

For a normal mode wave disturbance of the form 

p = Re {\Jj(y,z) exp i k(x - ct)} , (6) 

where ^p(y,z) is the perturbed pressure, k is the zonal wave number 
(taken as real) , and c is the phase speed (taken as complex , c = c + 
ic.) ; equations (2), (4), and (5) become, respectively, ^ 

(u - c) + ipyy - = 0 , (7) 

(u - c) - u i|j = 0 , on z = 0, 1, (8) 

z z 

i|;=0,ony = 0,l, (9) 

where the suffixes denote differentiation. 

To solve the problem we take u as a function of z only and put 
8 = 0. Equation (7) then takes the form 

(u - c) - (Fu^)^ = 0 . (10) 

Equation (10) is separable. Assuming a solution of the form ip(y ,z) = 

X (y) 4>(z) and using the boundary conditions (9), we obtain 


’jj(y ,z) 


(p(z) sin miry 


( 11 ) 


0 


( 12 ) 


(u - c) {(F())^)^ - K2<j,| - (Fu^)^ 4 , = 

2 2 2 2 

where K =k + niTT,m = l, 2 The boundary conditions (8) 

become 


(u - c) (j) - u (|) = 0 , on z = 0, 1 . (13) 

z z 


We are concerned with the solution of equation (12) for variable 
gravity. For a general power law, equation (1) can be written in the 
form 


g 



(14) 


where g^ and are values at the origin of the coordinates and n is a 

constant. Using equation (14) and assuming that all other quantities in 
F are independent of z, we obtain 


F(z) = 



(15) 


9 2 2 2 1 /9 

where F = f L /N H and N = (vg 3T/9z) and where F and N 
o o o 00 

are independent of z. From the thermal wind relationship, we have for 
the dimensional vertical shear of the basic state 


U^(z) 


gY 


dT 

3y 


Using equation (14) and again assuming that all other quantities are 
independent of z, we obtain 


Uz(z) = 



(16) 


5 



where = g^yf 


-1 


8T/8y and c is a characteristic value of U (z) used 
y z 

to nondimensionalize U^(z). and are independent of z. Let us 

consider the last term in equation (12). Note, using equations (15) and 

(16), that the combination Fu^ is independent of z, and hence this term 

z 

vanishes identically. Equation (12) then reduces to 


(F<f.^)^ - = 


(17) 


Expanding equations (14), (15), and (16) to first-order in y, 
we obtain 


g = go " 

yz) 

(18) 

F = F (1 + 
0 

yz) 

(19) 

-1 



6 u = (1 

z 

- yz) 

(20) 


where y = nH/z and 6 
to z we obtain 



Integrating equation (20) from z = 0 


6 




( 21 ) 


Equations (18) through (21) approximate the more exact forms previously 
given when y is small. In the next section we shall solve this approxi- 
mate variable g problem. 

The solution of the problem for constant g, (i.e. y =0) was 
first given by Eady (1949). A succinct statement of Eady's problem has 
been given by McIntyre (1970). Eady's problem is formulated from the 
above equations as 

(p - K^<i> = 0 , ( 22 ) 

ozz ^o 


(u-c)d)-U(}) =0 , on z = 0,l, (23) 

o o ^z oz^o 
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u = z , (24) 

o 

1/2 

where k = K/F . The solutions are 

o 

Cq= i±ia^{(a - coth a) (a - tanh a)} , (25) 

d) = K c cosh KZ - sinh k z , (26) 

^o o 


where a = k/ 2. The solutions represent either an ampUfying-^decaying 
pair of waves with = R^(c^) = 1/2, or a pair of neutri waves with 

> and < 1/2, according as k < or > , where 


= 2.3994 (k^ = coth k^) 


(27) 


There is a short wavelength cutoff to the instability at the critical 
neutral wavenumber . There is also a maximum growth rate for some 


III. SOLUTI ON B Y A PERTUR BATI ON METHOD 


In this section we shall solve the problem posed by equations 
(17), (18), (19), (13), and (21) using a perturbation method. The 
procedure followed is that described by McIntyre (1970) and it requires 

that p << 1. Figures 1, 2, and 3 show ^ ^ from 

equations (18), (19), and (21), respectively, plotted as functions of z 
for p = 0, 0.05, 0.1, 0.3 and 0.5. We take p to be a perturbation 
parameter and write 



= 1 + pF^(z) 

y 

(28) 

u(z) 

= z + pUj^(z) 

9 

(29) 
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Figure 1. Distribution of gravity, g/g = 1 - 
yz for the values of y as indicated. 


.2 


Figure 2. 
number , 




Distribution of the rotational Froude 

F/F = 1 + yz, for the values of 
o 

y as indicated. 




where 


F^(z) = z , (30) 

and 

Ui<z) = - I z^ . (31) 

The method is to expand (J)(z) and c about the Eady solutions (see 
Section II) in the following way 

(() = + y(|)^ + ■*■••• > (32) 

2 

c = c + yc. + y Cq + . . . . (33) 

O 1 z 


It turns out that this expansion, in integral powers of y , is 
appropriate only for the case k # where Kj^ is the transition wave- 

number from the neutral wave regime to the unstable wave regime. For 
the case k = the correction terms obtained from the expansions (32) 

and (33) diverge. This is because the Eady solutions possess branch 
points at K = Kj^. In view of this, two separate cases will be considered. 

For the first case k # k„ , the appropriate expansions are (32) and (33). 

^ 1/2 

For K expansion in powers of y has been found to be 

appropriate, namely 


1 /2 

<j) = (l)^ + y (f)j^ + y<})2+... > (34) 

1 /2 

c = c + y c. + yc„ + ... . (35) 

O J. z 


The substitution of equations (28) and (32) into equation (17) 
yields, to zero and first-order in y, respectively 
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0 


(36) 


d) - K d) 
^ozz 


♦izz - ^ ♦l 


Iz 


T? 


F ^oz F "^ozz 
o o 


(37) 


1/2 

where k = K/F^ . With the aid of equation (36), equation (37) becomes 


‘i’lzz “ ^ '*>] 


Iz 


, 12, 
F *oz ■ I- '' *< 


(38) 


In a similar fashion, the substitution of equations (29), (32), and (33) 
into the boundary conditions (13) yields, to zero and first-order, 
respectively 


(z - c ) (p - (p - 0 , on z = 0,l , 

o ^oz ^o 


(39) 


(z - c^) (j)^^ - (j)^ + (u^ - cp cf.^^ - u^^4>o = 0 , z = 0, 1 , 


(40) 


where c = 6 and c. = 6 . Combining these two equations, we 

obtain boll 


(z " c ) 
o 




0 , on z = 0, 1 . 

(41) 


Substituting from equations (30) and (31) into equations (38) and (41), 
we have finally 


Izz 


2 2 
-Kd), =-(}) -Kzd) 

^1 ^oz ^o 


(42) 
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. o^) - 


^ + 2g^ 

2(z - c ) 
o 


- z 


(f) =0,onz = 0,l 


(43) 


It can easily be shown that the general solution of the non- 
homogeneous equation (42) is 


= \ (k^z^ + 3c 

4 c 


2 

KZ + A) cosh KZ - (c IC z 

o 


+ 3z + B) sinh kz 


(44) 


which is the first-order correction to the eigenfunctions. With the aid of 
the boundary conditions (43) , the two constants A and B can be reduced 
to one arbitrary constant. Our primary interest, however, is the deter- 
mination of the first-order correction to the eigenvalue c. 

We multiply equation (42) by cf)^ and integrate from z = 0 to 

z = 1. The first term on the left is integrated by parts twice, and then, 
with the aid of (36), we obtain the result 




o^lz 


0 


- d> (p 
^oz^l 


/ ,2^ 1.2 

( Z(fl dz - jr 
Jr, O 2 


(45) 


0 


From the zero and first-order boundary conditions, equations (39) and 
(40) , we find 


6 ^(o) 
oz 


♦o(o) 


OZ 


1 - c 


cJ)l,(o) 


2 


(o) - 
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2(1 




[(2c^ + 


2c 


1 ) 


4>o(l) 


+ 2(1 


y tf>i(l)] 


Substituting these into equation (45), simplifying and solving for c. , we 
obtain ^ 


c 


1 




- / 

o 


1 


zd) (z) dz 
^o 


(46) 


With the aid of the Eady eigenfunctions (26) , this expression for Cj^ can 

be expressed in terms of hyperbolic functions. A neat form can be 
derived by eliminating the hyperbolic functions using identities derivable 
from equation (25). After some substantial and judicious algebraic 
manipulation, it can be shown that 


(1 - 


c. = 


~ . 2,^2 2 

Oq) (c^ k - 1) 

4(1 - 2c^) 


o 

4 


(47) 


Equation (47) involves only the unperturbed eigenvalue c^. We 
note that if c^ is real, c^ is real. If c^ is complex, the real and 
imaginary parts of c^, are given by 


^ir 

1 o 2 , 

o ( R 3 cx ) 9 

(48) 

16a 



1 rr 2 2 _2. 


°li 

~ 16^ + R)-a - R] 

(49) 

2 

where R = 

- (a - coth a) (a - tanh a) and a = k/2. 

For c real, c, is 

given by 


o 1 
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Ir 


li 


16a 
0 , 


1 /o2 ^ O 2. . 1 2 „2. ,2 c.2^2, 

(S + 3a ) — ^gojg t(oi S ) (a “* S ) J $ 


(50) 


2 2 

where S = - R . Thus, for a < coth a, (k < Kj^) we have finally. 


I + (R^ - 3a^) 

^ IGa'^ 


(51) 


_ . . / R ^ u r/ 2 „2. ,2 „2.2, \ 

T “ I 2a 16oR ^ ^ R ) ] | 


(52) 


and for a > coth a, (k > Kj^) > 


c 

r 


c. 

1 


6 

2 

0 


(S^ + 3a^) ± 6 

16a 


S u 
2a ^ 16aS 


r. 2 _2, ,2 o2,2J 

[(a - S ) - (a - S ) ]| , 


(53) 

(54) 


The quantity 8 (= ^) depends on a characteristic value for 

U (z). We shall first plot the results for C = C > and hence for 6=1. 
z ^ p o 

This means that for the characteristic value of c we have chosen the 


single representative value at z = 0. Figure 4 shows the phase speeds 
calculated from equations (51) and (53) and plotted as c^ versus k for 

y =0.1, 0.3 and 0.5. The Eady result (see Section II) corresponding 
to y = 0 is also shown. Figure 5 shows the growth rates calculated 
from equations (52) and (54) and plotted as kc. versus k for y = 0, 

0.1, 0.3 and 0.5. Note that the values diverge as the wave number 
approaches the Eady critical value k = = 2.3994. 


As was discussed earlier in this section, the expansions (34) and 
(35) must be used at k = Kj^. Following the general procedure used 

above for k # it can be shown that for k = 


14 



Figure 4. Variation of the phase speed, c^, for 6 =1, with the wavenumber k 

for the values of y as indicated. 


2 


K 


3 


4 


normalized growth rate, kc^, for 6=1, with the wavenumber k 
for the values of y as indicated. 



( 55 ) 




K (K 2 - 4)^/2 

16 


Thus, we have 


°r 2 


K (<2_4)1/2 

^ 1/2 ^ N '• N ^ 

c. = ±u S —16 


(56) 


(57) 


(58) 


These results are shown as dots in Figures 4 and 5 for 6=1 and for 
the same values of y as the curves already plotted. 


The plots in Figures 4 and 5 are based on representative values 
of F and u^ at z =0. This means that they show the effects of both an 

average lower value for g as well as the variation of g, whereas plots 
based ion characteristic average values of F and u^ would show only the 

consequences of the variation of g. Let us choose linear averages over 

z for F and u , defined by 
z 


F 


and 


f F(z) dz 
0 


6 

y 



u (z) dz 
z 


Substitution from equations (19) and (20), respectively, and performing 
the simple integrations , we obtain 


F = 



implying that k 


K(1 



(59) 


17 


6 


(60) 



Figures 6 and 7 present plots of the phase speed and growth 
rates based on equations (51) and (53), and (52) and (54), respectively, 
and using equations (59) and (60). The dots are based on equations 
(57) and (58) and using (59) and (60). 

Since only the first-order corrections are calculated, the results 
are accurate only for small y . The results should be reasonably valid 
for y = 0.1, but for the larger values of y moderate errors are certainly 
present. 


IV. DISCUSSION 


In this paper we formulated an approximate, quasi- geostrop hie , 
baroclinic stability problem for a body force with a height dependence. 
The primary motivation for this work was to determine the effect of a 
variable dielectric body force, analogous to gravity, on baroclinic 
instability for the design of a spherical, atmospheric, synoptic-scale 
model experiment for Spacelab. In the formulation, curvature and 
horizontal shear of the basic state were omitted and the vertical and 
horizontal temperature gradients of the basic state were taken as con- 
stants independent of z. A consequence of the body force variation and 
the above assumptions is that the transverse gradient of the potential 
vorticity of the basic state vanishes. This stability problem was then 
solved using a perturbation method. 

Figure 5 shows the well known results due to Eady (y =0), that 

for K = Kxt = 2.3994, unstable waves exist, and for k > k,, neutral waves 
N N 

exist (see Section II). The growth rates for the non-zero values of y 

are qualitatively similar to Eady's results and converge to Eady’s result 

as y -> 0. Note that, because of the divergent behavior at » we 

cannot accurately locate the critical value of k for each value of y . 
However, by extrapolation using the point values, we can see that the 
critical values move to larger values of k as y increases. Figure 5 shows 
that in general the growth rates decrease as y increases. This is not 
surprising since a dimensional statement shows the growth rate depends 
linearly on g, and g decreases as y increases. Note the cross-over 
point to the right of which the opposite is true. 

Figure 4 shows that the phase speed of the unstable mode 
decreases as y increases and is only weakly dependent on k. Using 

equation (21) we obtain for the average value of 6 ^u, 6 ^u = 1/2 - y/6. 
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For 6=1 and y =0.1, 0.3 and 0.5, u = 0. 483, 0.450 and 0.417, 
respectively. The phase speeds are in good agreement with these values. 
Thus, like the Eady result, the unstable mode is simply advected at the 
rate of the mean flow. The levels at which the actual flow corresponds 
to the phase speed of the unstable wave are known as the steering levels, 
and are found by solving 


z 

s 



2 

z 

s 



In the stable region there are two neutral modes, one confined near the 
upper boundary and one near the lower boundary (Bretherton, 1966; 

Hyun and Fowlis, 1979). Each of these neutral modes also propagates 
with a phase speed related to a mean flow in its environment. The phase 
speed of the lower mode does not vary much as a function of y since the 
basic flow does not change much in the lower region (Fig. 3). However, 
in the upper region there is a substantial change as a function of y in 
the flow, and this is reflected in the phase speed. 

Figures 6 and 7 show the phase speeds and growth rates plojtted 
in terms of averaged values of the basic state. For small values of k 
the growth rates are now almost independent of y but for values of x 
near the maximum growth rate and near the critical values, the growth 
rates are larger. 

We now proceed to ascertain what information the calculation has 
for the design of the Spacelab experiment. For this experiment, n = 5 
[see equation (1)1 and the gap-to-mean-radius ratio, H/z^, will be 

approximately 0.1. Thus y = nH/z^ - 0.5. Our first-order perturbation 

calculation cannot be expected to give an accurate answer for such a 
large value of y , but an approximate result and the general behavior of 
the instability as g varies can be expected. [We should also remember 
that the binomial approximation was used to formulate the equations (18), 
(19) and (20) and the accuracy here depends again on y being small.] 

The plots show that the results are qualitatively similar to Eady^s results. 
No new phenomena have appeared. The growth rates based on averaged 
values are somewhat larger than for the Eady problem, and the range of 
instability is somewhat increased. Phase speeds remain close to mean 
flow speeds. Thus, the results suggest that for a body force variation 
of 50 percent, the stability characteristics are not much changed. Future 
design studies should treat this problem more accurately. 

Although the effect of variable gravity on baroclinic instability, 
per se, is not of much value in understanding atmospheric and oceanic 
flows and laboratory geophysical fluid flow experiments, the solutions we 
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have obtained can be relevant to these flows. The variation in g which 
leads to a variation in the Brunt-Vaisala frequency, N, and in turn to 
a variation in the rotational Froude number, F, could formally have been 
due to a variation in the vertical temperature gradient , 8 T / 8 z , with g 
constant (see the definitions of N and F in Section II). Thus, if we 
deal with a flow for which N is a function of z and for which the vertical 
shear is also a function of z, such that the gradient of the potential 
vorticity, q^., vanishes, then our solutions are relevant. The ocean and 

the laboratory cylindrical annulus (see Section I) are systems in which 
the static stability undergoes large variation with height. 
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